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Introduction 

The modified Langley 8 -Foot High- Temperature Tunnel ( 8 ’ HTT) will utilize pressurized gaseous nitrogen 
(GN 2 ) at about 2000 psig to transfer liquid oxygen (LOX) from a run tank to the tunnel combustor. Since 
the contamination of the LOX by pressurized GN 2 may prove to be an important operational constraint, it is 
desirable to have reliable data concerning the penetration of nitrogen into LOX during pressurization. Such 
data, unfortunately, are scanty or nonexistent. 

Contamination of the LOX takes place through the following process. When GN 2 is first admitted to the run 
tank, it cools down rapidly to near the LOX temperature (« 90.2 K). When the pressure reaches approximately 
52 psia, a layer of liquid nitrogen (LN 2 ) builds up at the LOX interface. The LN 2 penetrates into the LOX 
primarily by simple molecular diffusion. The purpose of this study is to determine experimentally the extent 
of the penetration of LN 2 into LOX under given operating conditions. 

The principle of the measurement is based on a technique known as “differential flash vaporization” (ref. 1), 
whereby a two-component liquid is permitted to evaporate and the vapor formed is removed from the remaining 
liquid. An analysis of the vapor yields sufficient information to determine the initial composition of the liquid, 
from which the binary diffusivity can be obtained. The adaptation of this technique to the present experiment, 
which was conducted in one of the 4-gal run tanks of the Langley 7-Inch High- Temperature Tunnel (the pilot 
tunnel for the 8 ’ HTT), is illustrated in figure 1. After the run tank is filled with LOX to a specified level, the 
test consists of three basic steps: 

1. Pressurization — the run tank is pressurized with GN 2 for a specified holding time 

2. Blowdown — the pressurized GN 2 is vented to 0 psig 

3. Analysis — the vapor from the evaporating liquid is analyzed in an oxygen analyzer, and a time history of 
the O 2 content is recorded for subsequent data processing. 



(a) Pressurization 



(b) Blowdown 



Figure 1. Principle of the differential flash vaporization technique. 

This paper is organized into three basic parts: theory, experiment, and results. The theory includes a 
physical model of the LOX contamination process and a derivation of the equations used to evaluate the 
diffusivity from the data. The experimental sections contain first an account of a preliminary experiment, 
using a nuclear monitor, to study the buildup of the LN 2 layer under pressurization, and then a description 
of the experimental procedure used in the differential flash vaporization experiment. The results include the 
evaluation of the diffusivity, error analysis, application to operating conditions in the 8 ’ HTT, and comparison 
with prior measurements at lower temperatures. 

The authors express their thanks to Andrew P. Seamons, Governor’s School for the Gifted, 1987, for 
assistance in making the measurements. 

Symbols 

A surface area of liquid, cm 2 

b \ , 62 > &3 parameters to fit pressure time history 

D binary diffusivity of LOX-LN 2 mixtures, cm 2 /sec 
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preexponential factor, cm 2 /sec 

diffusion length (penetration depth), cm 

total number of moles of liquid, mol 

total number of moles of vapor in ullage volume, mol 

total pressure in the run tank, psig 

reference pressure for flow rate calibration, atm 

pressure at the flowmeter, atm 

evaporation rate, mol/min 

universal gas constant, cm 3 -atm/mol-K 

reference flow rate, cm 3 /min 

temperature of vapor in ullage volume, K 

temperature of the liquid, K 

reference temperature for flow rate calibration, K 

temperature at the flowmeter, K 

time, min 

holding time, min 

time at inflection point, min 

ullage volume, cm 3 

molar volume of liquid, cm 3 /mol 

molar volume of LN2, cm 3 /mol 

molar volume of LOX, cm 3 /mol 

mole fraction of N2 in the vapor 

mole fraction of N2 in the liquid 

value of x nl fi me t — 0 

location along the run tank, in. 

distance from interface, cm 

interface location at any time, cm 

interface location at time t = 0 , cm 

activation enthalpy, kcal/mol 

coefficient of volume expansion, K _1 

liquid density, g/cm 3 

time constant to flush out ullage volume, min 
characteristic evaporation time, min 
value of t\ at unit pressure, min 

indicates inflection point 
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Superscript: 

' (prime) indicates quantity for holding time of 1 min 

Mathematical functions: 

2 

erf(x) error function, (2/^n) Jq e~ u du 

erfc(x) complementary error function, 1 - erf(x) 

erf^^x) inverse error function 

Note on units: Quantities read directly from tunnel instrumentation retain the units of the instrumentation, 
usually English units. Quantities used in the evaluation of the data are expressed in metric units. 

Theoretical Background 

Physical Model of the LOX Contamination and Monitoring Processes 

The model is described with the aid of figure 2, showing the N2 profile during pressurization, blowdown, 
and analysis. Prior to pressurization the run tank is filled with LOX to a specified level. After pressurization a 
layer of LN2 builds up at the LOX interface. The concentration profile of LN2 immediately after pressurization 
is indicated by curve A in figure 2(a). The origin 2 = 0 is chosen as the LOX-LN2 interface location. Initially 
the run tank contains pure LOX for 2 > 0 (toward the bottom) and pure LN2 for 2 < 0 (toward the top). 



(c) During analysis. 



Figure 2. Distribution of LN 2 in the run tank during the three steps of the experiment. 


The large concentration gradient across the interface results in mass transport by molecular diffusion. After a 
given holding time the LN2 concentration profile is given by curve B, indicating that the LN2 has penetrated 
into the region formerly occupied by the pure LOX. The quantity of LN2 passing into this region increases 
with increasing holding time. Curve B, representing a typical LN2 concentration profile at the end of the 
pressurization step, is described mathematically as a complementary error function (ref. 1): 

^=GMvk) = G) 

where D is the diffusivity and tf t the holding time. 

During blowdown the falling pressure quickly reaches a level that can no longer sustain N2 in the liquid 
state, and the N2-rich portion of the liquid immediately evaporates. When the pressure finally reaches 0 psig, 


1 — erf 


(2 JTnJ_ 


( 1 ) 
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the liquid has evaporated to some point zq below the original interface location as the result of N 2 penetration 
into the LOX. The concentration profile immediately at the end of blowdown is illustrated in figure 2(b). At 
this time (defined as t = 0) vapor from the evaporating liquid is routed to the oxygen analyzer. 

During the analysis step the interface location Z{ shifts steadily from its initial location zq toward the 
bottom of the run tank (increasing z). At any time t > 0 the composition of the vapor is representative of that 
of the interface location from which it has evaporated. The steadily progressing shift in interface location is 
illustrated in figure 2(c). 


Simplifying Assumptions 


A mathematical description of the events depicted in figures 1 and 2 is intractably complex but can be 

made feasible with the aid of several simplifying assumptions: 

1. During pressurization there is no thermal (Soret effect) or convective mixing of the liquids. The temperature 
difference between the LOX and the LN 2 at the interface is too small to provide a significant driving force 
for thermal diffusion. Convective mixing ensues from an instability developing between two opposing effects. 
The natural buoyancy due to the lower density of LN 2 compared with LOX, 0.7912 vs. 1.167 g/cm 3 (under 
pressurization), is a stabilizing effect. The lateral temperature gradient, due to the temperature difference 
between the wall of the run tank (77.4 K) and the core of the liquid (90.2 K), is destabilizing in that it 
causes thermal buoyancy which tends to initiate a circulating axial flow. If one adopts the criterion that 
the thermal buoyancy must exceed the natural buoyancy for convective mass transfer to occur, then the 
thermal Grashof number must exceed the concentrative Grashof number (see ref. 1, p. 66). The ratio of 
these numbers is p^K ATl/ A pi, which in the present case has the value, in cgs units, 


1.167 x 0.0038 x 12.8 
1.167 - 0.7912 


= 0.15 < 1 


indicating a high degree of stability against convective mass transfer. 

2. There is no net molecular transport across the interface; that is, for every N 2 molecule that crosses from 
above, a corresponding O 2 molecule crosses from below. The practical significance of this assumption is 
that the interface location shown in figure 2(a) does not change during pressurization. This manner of 
transport is known as “equimolar counter diffusion,” which is characterized by a single coefficient known 
as the “binary diffusivity.” 

3. The concentration profile of the liquid remains frozen after blowdown. The LN 2 remaining after blowdown 
in figures 2(b) and 2(c) follows the profile at the end of the pressurization interval, curve B in figure 2(a). 
The justification for this assumption is that the concentration gradient is too small, except at the rapidly 
evaporating leading edge, to permit significant diffusion after blowdown. 

4. Evaporation of the liquid proceeds at a constant rate R. Actually, as will be shown later, this condition 
can be enforced since the evaporation rate can be adjusted by means of a hand valve. Consequently, the 
interface location Z{ advances linearly with time: 


Zi = zq + Rv L t/A (2) 

where R is the evaporation rate, vi the molar volume, and A the cross-sectional area of the liquid. 

5. The steady-state vapor flow rate is constant and equal to the evaporation rate. In other words, after 
blowdown, pressure does not build up in the run tank. As will be seen, this assumption is not quite valid, 
but a method for correcting for a slow pressure buildup is described in appendix B. In any case, a breakdown 
of this assumption has minimal impact upon the evaluation of the diffusivity. 

6. The vapor in the ullage volume, i.e., the volume of the run tank above the liquid plus that of the system 
piping, is homogeneous. Vapor from the evaporating liquid mixes very rapidly with that already present in 
the ullage volume. The O 2 analyzer samples the homogeneous vapor mixture. 

7. The ullage volume is constant. Because most of the diffused N 2 is contained in a small volume near the 
interface, the essential information concerning LN 2 concentration is gathered from the initial small fraction 
of the evaporating liquid after blowdown. 
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Analysis 

In figure 1(c) let ny be the total number of moles of vapor in the ullage volume and the total number 
of moles of liquid. Conservation of molar content requires, as the liquid evaporates, that 


Similarly, conservation of moles of N 2 requires that 

J t ^ x N n v) = -~jj.(xNL n L) 

^ + ( 4 ) 

where xj v and x^l are the mole fractions of N 2 in the vapor and liquid, respectively. At this point the 
analysis takes a different course from that taken in conventional flash vaporization. In the latter the liquid is 
a homogeneous mixture of two components. As the liquid evaporates, the mole fraction of the more volatile 
component in the vapor changes very slowly. Thus dx^/dt is assumed small, and the second term in equation (4) 
is neglected. 

In the present experiment the situation is considerably different. The liquid is not a homogeneous 
mixture; rather, x^i is strongly spatially dependent and x^y changes rapidly with time. As a consequence of 
assumption (3) (frozen LN 2 profile), on the other hand, x^i does not change with time: 

< 5 > 

and the last term in equation (4) vanishes. Substitution of equations (1), (3), and (5) into (4) leads to the 
following transport equation at the interface z = zy. 

= [ l - ert (_^)] ( 6 ) 


Upon defining, for convenience, 


h = 2 y/Dth 

ny 

Ti = -r 

Ah 

r 2 = 77 — 


and using equation (2), one can rewrite equation (6) in more simplified form: 

^v + 5v = a?i= 1 f 1 _ erf ('2L + i')l (io) 

dt t\ T\ 2t\ \ h t 2 J \ 

Each of the convenience variables defined in equations (7)-(9) has a physical meaning: h is the diffusion length, 
or "penetration depth” corresponding to the holding time t,^: t\ is the time constant for flushing the GN 2 out 
of the ullage volume; and t 2 is a characteristic “evaporation time” representing the rate of advance of the 
interface. Noteworthy is the fact that, among the characteristic times, only r 2 contains h and thus information 
regarding the diffusion constant. 

The initial condition is expressed in terms of the mole fraction of N 2 in the liquid at the end of blowdown, 
in other words, at the time the analysis begins. Using equations (1) and (2), one expresses this condition as 


= x N d °) = \ 1 - erf 
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This is indicated in figure 2(b). Details of the solution to equation (10), together with initial condition (11), 
are given in appendix A. The solution is 


[l — erf (f + -)1 

1 

+ x exp 

L \ h T2 / J 

2 


L + Z ° T 2 

T\ hr\ 
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e^t (4 + r"^)■' !^,| 


(?-£)] (“> 


In sum, equation (12) represents the time dependence of the mole fraction of N 2 in the ullage volume — in other 
words, the mole fraction that is detected in the oxygen analyzer. 

Preliminary Investigation of the LOX-LN2 Interface Profile 

A test to gather information concerning the structure of the LOX-LN 2 interface profile was conducted using 
a nuclear monitor developed by Singh, Davis, and Mall (ref. 2). This monitor is capable of identifying the 
medium — LOX, LN 2 , or vapor — at any location within the run tank. 

The test setup is illustrated in figure 3. A Cs 137 source emits a collimated beam of 662 keV 7 -radiation to 
a Nal detector on the opposite side of the run tank. The detector counts the number of emissions within a 
given time interval. The beam is attenuated by the steel wall of the run tank as well as by the medium within 
the wall. Despite a rather high attenuation in the former, the detection system is still sensitive enough to 
distinguish among attenuation differences in LOX, LN 2 , and vapor. Table I shows the number of 15-sec counts 
with each of these three media contained in the space within the wall, together with their standard deviations. 
It is clear that the count yields an unambiguous identification of the medium. The source and detector are 
mounted on a rigid frame, which slides along a track. The frame can be locked at any location Z, as measured 
by a scale along the tank from 0 to 100 in. The tank extends an additional 28 in. beyond either end, regions 
which are not accessible to the monitor. 



Figure 3. Test setup for the preliminary investigation of the LOX-LN 2 interface profile. The 45° slant of the run tank is not shown. 


Table I. Number of 15-sec Counts of Cs 137 (662 keV) Gamma Rays Passing Through Various Media 


Medium 

Number of counts 

LOX 

6700 ± 161 

ln 2 

8300 ± 161 

Vapor 

10 500 ± 95 


In the interface profile tests the nuclear monitor was fixed at a location Z = 50 in., halfway along the run 
tank. The run tank was precooled with LN 2 , emptied, and filled with LOX to the 50-in. level, as determined by 
a change in count on the monitor. After GN 2 was admitted to the run tank to a predetermined pressure, the 
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monitor was moved to identify the medium at selected locations starting at Z = 30 in. A typical result is shown 
in figure 4 for a GN 2 pressure of P — 750 psig. At each location a set of three 15-sec counts was recorded, and 
the monitor was then moved to the next location. After the measurements at Z = 100 in. (top of the track), 
the monitor was moved back to Z — 55 in. to confirm the previous readings. The count indicates the presence 
of LOX between 30 and 50 in., a transition region from 50 to 55 in., and a layer of LN 2 between 55 and 70 in. 
At locations of 85 and 100 in. the count is too high for LN 2 but too low for vapor. The interpretation of 
these values is based on the fact that the run tank becomes warmer near the top. At locations above 70 in. 
the temperature apparently exceeds the critical temperature, and N 2 exists in the supercritical state with 
decreasing density toward the top of the run tank. 



Location, in. 


Figure 4. Interface profile at a GN 2 pressure of P = 750 psig. 

The test was repeated for GN 2 pressures of 250, 500, 1200, and 1800 psig. The corresponding layers of LN 2 
above the interface measured 3, 10, 30, and 40 in., respectively. For the present discussion, N 2 at temperatures 
near 90 K is considered a liquid even at pressures beyond the critical pressure of 492 psia. In conclusion, under 
pressurization a layer of LN 2 builds up at the LOX interface to a level where the temperature is too high 
to support the liquid. Since density at the interface conforms to that of the liquid rather than the vapor, a 
consequence of the LN 2 buildup at the interface is enhanced mass transport through molecular diffusion. 

Experimental Procedure 

The LOX-LN 2 contamination experiments were carried out at the Langley 7-Inch High- Temperature Tunnel, 
the pilot tunnel for the Langley 8-Foot High- Temperature Tunnel. The experimental setup is shown in figure 5. 
The piping in the system, with the exception of coils 1 and 2, was oxygen-cleaned and inspected. The tests 
consisted of four major steps: (1) setup, (2) precooling of the run tank, (3) pressurization and blowdown, and 
(4) monitoring. 

During the setup step of each test the monitoring instrument, a Beckman Model 755 oxygen analyzer 
(0-100%), was calibrated, and the nuclear monitor was installed. Calibration procedures required the use of 
the 3-way valve D shown in figure 5. Pure bottled GN 2 (99.995%) was used to calibrate the zero point, and 
bottled air was used to set the 20.9% O 2 point. The Beckman analyzer was checked on one occasion against 
uncontaminated evaporating LOX in the 4-gal run tank and yielded a reading of 99.992%. The nuclear monitor 
was installed at the 46-in. mark on the scale outside the run tank to give a liquid level indication for oxygen 
filling. 

Since the thermal insulation of the run tank is not of Dewar quality, two methods were implemented to 
slow the evaporation of the LOX-LN 2 mixture. First the run tank was precooled by filling it completely with 
LN 2 . After an ample cooling period, the LN 2 was discharged to the atmosphere. Second, LN 2 was injected 
through two copper tubes that coiled around the outside of the run tank as seen in the figure. The addition 
of the cooling coils decreased the evaporation rate by a factor of 5 compared with earlier experiments. At the 
end of the precooling period LOX was transferred immediately into the run tank and filled to the 46-in. level 
as measured by the nuclear monitor. Normally filling continued for 1 additional minute after this point was 
reached to ensure that the liquid level stabilized. 
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Figure 5. Test setup for the differential flash vaporization experiments. 

The next step, pressurization and blowdown, began with pressurizing the top of the run tank with GN 2 . 
Various GN 2 pressures and holding times, listed in table II, were used to study how these parameters affected 
the N 2 penetration into the LOX. After the GN 2 was held for the predetermined time, the pressure was vented 
to prepare for monitoring of the evaporating LOX-LN 2 through the Beckman analyzer. 

Table II. Summary of Experimental Parameters 


Run 

N 2 pressure, psig 

Holding time f/,, min 

Comment 

1 

1700 

15 

Supercritical for both O 2 and N 2 

2 

1700 

15 

Supercritical for both O 2 and N 2 

3 

1700 

15 

Supercritical for both O 2 and N 2 

4 

1700 

15 

Supercritical for both O 2 and N 2 

5 

1700 

15 

Supercritical for both O 2 and N 2 

6 


15 

Supercritical for both O 2 and N 2 

7 

500 

15 

Supercritical for N 2 but not for O 2 

8 

1700 

1 

Supercritical for both O 2 and N 2 

9 

30 

15 

Normal 


For the monitoring step the pressure in the run tank was vented to 0 psig before venting through the 
Beckman analyzer to ensure that the residual GN 2 was flushed out of the monitor line. After valve 3409X was 
opened the flow from the run tank proceeded through two parallel paths. The first path passed the bulk of the 
vapor through flowmeter 1 at 30 000 cm 3 /min, as controlled by valve A, and was then vented to the atmosphere. 
The second path passed a sample through the Beckman analyzer and flowmeter 2 at 250 cm 3 /min, as controlled 
by valve B, and then was also vented to the atmosphere. The flowmeter readings were continuously adjusted to 
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Figure 6. Time history of N 2 content in the vapor monitored by the Beckman O 2 analyzer for a typical run (Run 2). The GN 2 
pressure was P = 1700 psig for a holding time of t h = 15 min. Three regions are identified: I (0-1 min), adjustment time 
of the Beckman monitor; II (1-4 min), evaporation of the N 2 -rich liquid; and III (> 4 min), flushing out of GN 2 from the 
ullage volume. 

these fixed levels on each test run. Valve C was left open during each run and served only to prevent backflow 
from the atmosphere when the system was not in use. 

For the remainder of the run, data were recorded that included the temperature of the cooling coils as 
measured by a thermocouple, the percent O 2 as measured on the Beckman analyzer, the pressure in the 
ullage volume, and the elapsed time. With the exception of the pressure, which was recorded manually, these 
measurements were recorded in 10-sec intervals by means of an Omega analog-to-digital converter interfaced 
to an Apple He computer system. The experiments ran until the Beckman analyzer reading neared 100% O 2 
for a period of time. At the end of the run the Beckman analyzer was recalibrated. 

Experimental Results 

The first five runs listed in table II were taken at the maximum available GN 2 pressure of 1700 psig for a 
holding time of 15 min. It is known from experience that within this time conditions remain reasonably stable 
inside the run tank. The next two runs were conducted at reduced pressures, at the same holding time, to 
test for possible pressure dependence of the diffusivity. The next run was conducted at a shortened holding 
time to test for the known time dependence of penetration by diffusion. A final run was conducted at a GN 2 
pressure of 30 psig, which is too low to establish a layer of LN 2 . Since diffusion is assumed to be negligible in 
the absence of a liquid layer, this run permitted the response time of the oxygen analyzer to be determined. 

Experimental Time History of the N2 Mole Fraction 

The time history of the N 2 content of the vapor passing through the oxygen analyzer, shown in figure 6 for 
a typical test (Run 2), is marked by three distinct regions: I, a sharp drop from 100% to a little under 40% 
in the first minute; II, a small crest in the interval 1-4 min; and III, a slow decay for times exceeding 4 min. 
These are shown more clearly in figure 7, which gives an expanded view of the first 10 min. The transition 
from one region to another indicates a gradual change in the physical process governing the time history; the 
lines of demarcation serve only to illustrate the approximate times when one process yields to another. 
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Figure 7. Expanded view of the first 10 min of figure 6. The point where the derivative has an extreme value is indicated by 

“M.” 

The sharp drop in region I represents the response time of the O 2 analyzer, as opposed to the transit time 
delay, for the transit time of the vapor from the run tank to the analyzer is estimated to be a rapid 0.2 sec at 
a flow rate of 30000 cm 3 /min. At the start of the analysis step (see fig. 1), the O 2 analyzer was full of pure 
N 2 , left over from the earlier calibration. After the sample inlet valve was opened, the analyzer required a 
finite time to adjust to the prevailing value of the test gas. According to the instruction manual (ref. 3), the 
instrument requires 20 sec to reach 90% of the steady-state reading, following a step change in concentration. 
This figure corresponds to a time constant of 8.7 sec. In the present experiment the adjustment interval was 
chosen to be 5 datum sample times, or 60 sec, for this is the time it takes for the analyzer reading to reach 
[1 — exp(— 60/8.7)] = 99.9% of its steady-state value. The third significant figure is within the resolution of the 
analyzer. 

The above interpretation of region I was tested in a specially designed experiment (Run 9), in which the 
GN 2 pressure was made too low (30 psig) to permit the buildup of a liquid layer at the interface. Consequently 
there was negligible penetration of N 2 into the LOX and, after blowdown, the evaporating liquid consisted 
of pure O 2 . The resulting time history is shown in figure 8. Here a very rapid but finite decay of the N 2 
mole fraction is evident. At the fifth datum point the N 2 mole fraction measures 0.7% (99.3% O 2 ), in good 
agreement with the value computed above. Thus this test confirms the existence of an adjustment interval and, 
as a result, the first five data points have been excluded from the evaluation of the remaining experimental 
runs. 

The small crest appearing in region II is attributed to the evaporation of the N 2 -rich portion of the liquid. 
In the example of figures 6 and 7 this occurs in the interval 1-4 min, when the mole fraction of N 2 drops from 
about 40 to 35%. At the end of this region nearly all N 2 is out of the liquid. 

During the slow decay shown in region III, residual GN 2 is being flushed out of the ullage volume. 
Continuous mixing of the residual GN 2 with vapor from the evaporating liquid, according to assumption (6), 
ensures that the GN 2 content never really goes to zero. 

A comparison of the experimental time history in regions II and III with the theory (eq. (12)) is beset with 
two major difficulties. First, the appearance of three adjustable parameters — zo/h, iq, and T 2 — in the nonlinear 
function given in equation (12) makes a least-squares fit to the data not only computationally intractable but 
of questionable reliability, for the estimate of one parameter is sensitive to small errors in the estimates of the 
other two. Second, the parameter t\, which is proportional to the number of moles in the ullage volume, is 
thus a function of temperature, pressure, and volume. The fact that these parameters are not constant over 
the full duration of the decay causes assumptions (5) and (7) to be violated. 
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Figure 8. Time history of N 2 content in the vapor monitored by the Beckman O 2 analyzer for a low pressurization run (Run 9). 
The GN 2 pressure was P = 30 psig, too low for condensation to the liquid state, for a holding time of t/ t = 15 min. 

A simplified method of evaluating the adjustable parameters has been devised to avoid these difficulties. 
First, the parameter zq/H is determined from the initial condition, equation (11), where time t — 0 is defined 
as the time of the fifth datum point. In other words, the theory is forced to fit this point. 

Second, the parameter t\ is eliminated by the following technique. Differentiation of equation (10) yields 

d 2 XN 1 dxpf _ 1 (t/ T7+ z n /h) 2 / 10 ') 

dt 2 + T\ dt ^TJT 2 ^ ^ 

Now, at the inflection point M (see fig. 7), t — t^, 


'd 2 x 
dt 2 


N 


= 0 


M 


and dxj^/dt assumes an extreme value (dx jy / dt) m . Equation (13) becomes independent of t\\ 


f dx N \ 

\ dt ) M 


e -( J Af/r2+W ft ) 2 


\f^ T 2 


(14) 


Since (dx^/dt)M and t^ are known from the data, equation (14) yields T 2 - The values of T 2 derived from 
equation (14) are listed in table III. 

Table III. Summary of Results for Runs Used to Evaluate the Diffusivity 
[N 2 pressure = 1700 psig, holding time = 15 min] 


Run 

x m 

z 0 /h 

V£, cm 3 /mol 

T 2 , min 

1 

0.171 

■Mil 


0.91 

2 

0.394 



0.79 

3 

0.260 



0.92 

4 

0.299 

0.3915 

29.00 

0.74 

5 

0.181 

0.6443 

28.64 

0.88 

Mean ± standard deviation 

28.90 ± 0.269 

0.848 ± 0.0708 
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Finally ri, as defined in equation (8), is directly proportional to ny, which is given by the ideal gas law: 


ny 


PV 

RgT 


(15) 


In the early period of evaporation V and T are assumed to remain constant, but P builds up very slowly in the 
ullage volume. The measured time dependence of the pressure for each of the runs can be fitted to an equation 
of the form 


P — 14.696 b\ + b 2 t — (&i — l)e b3t 


(16) 


In Run 2, for example, a fit to the pressure data yields b\ = 6.86, b 2 = 0.089827, and 63 = 0.09. If t is expressed 
in minutes, then equation (16) gives P in psia. Substitution of equations (15) and (16) into equation (8) leads 
to an expression for t\ of the form 


n = T R 


61 + b 2 t - ( bi - \)e~ b3t 


(17) 


where tr is the value of t\ at time t = 0. Upon substituting equation (17) for t\, and inserting the values of 
za/h and t 2 obtained from the above procedures, one finds that a fit of equation (12) to the data of figure 6 
yields tr = 7.8 min. The resulting theoretical expression is compared with experiment in figure 9. The 
fit is considered excellent up to the time t = 20 min, after which time the experimental values veer away 
from the theory. The reason for the discrepancy is believed to be a breakdown in the assumptions leading 
to equation (12), especially assumptions (5) and (7). Nevertheless, the excellent agreement during the initial 
20-inin period of the analysis step confirms the validity of the model, the assumptions, and the derived values 
of t\ and t 2 . 



Figure 9. Best fit of equation (12) (solid line) to the time history shown in figure 6. 
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Evaluation of the Diffusivity 

Equations (7) and (9) give the diffusivity in terms of the experimental parameters: 



( 18 ) 


In order to refer the diffusivity to specific conditions, only those five runs taken at a holding pressure of 
1700 psig and holding time of 15 min (Runs 1-5) were used in its final determination. The flow rate R , in 
mol/ min, is determined from 


cm 3 /min _ R 0 (T t P 0 jT 0 P x ) l / 2 ^ R 0 ( P 0 P \ \ 1 / 2 

cm 3 /mol R g T x IP x R g \ T a T\ ) { ’ 

where R 0 , P 0 , and T a are the reference flow rate, pressure, and temperature of the flowmeter during calibration, 
Pi and Xi the actual pressure and temperature during the test, and R g the universal gas constant. Substituting 
numerical values in equation (19) yields 


R = 


29 910 cm 3 /min 


lxl atm 2 


1/2 


82.053 cm 3 -atm/mol-K l 295.45 x 297.15 K 2 


= 1.230 mol/min 


( 20 ) 


The surface area of the liquid in the run tank is 

A = 38.3 x y/2 = 54.2 cm 2 (21) 

where the factor y/2 is needed because the run tanks are slanted at 45°, and the holding time is 

t h = 15 x 60 = 900 sec (22) 

The molar volume of the liquid mixture depends on that of each constituent weighted by molar content 
(assuming an ideal solution): 

V L = X NL V LN + (1 “ x NL) v LX = x Nl{^LN ~ V LX ) + V LX (23) 

At the start of region II, x^L = X;vo and at the end of region II, x^l ~ 0. The value of vi to be used in 
equation (18) is taken to be the mean of these two values: 

V L=\ [ x No( v LN - v Lx) + V LX + v Lx] 

= V LX + \ x N0( v LN - VLX) ( 24 ) 

According to reference 4, vix — 28.03 cm 3 /mol at T = 90.2 K. Since LN 2 is not a stable phase at this 
temperature at P = 1 atm, the value of LN 2 is assumed to be the same as that at T = 77.3 K, or 
V LN = 34.64 cm 3 /mol (ref. 5). Inserting these numerical values into equation (24) yields 

1 >i — 28.03 + 3. 30x^o (25) 

Thus vi differs somewhat from one run to another and is listed in table III. The mean values of vi and T 2 
over the five runs are used to evaluate the diffusivity. The numerical values given in equations (20)-(22) and 
table III yield for the diffusivity 

D = 8.6 x 10 -5 cm 2 /sec (26) 
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Estimate of the Error 


From equation (18) one can express the error in diffusivity D in terms of the errors in the experimental 
parameters: 

uni I X+. I I AD I I Am . I I I \KA\ 

(27) 


SD 



6 h 

±2 

6 R 

±2 

6 v l 

±2 

6 T 2 

±2 

SA 

D 


t h 


R 


vl 


72 


A 


where v i and T 2 are mean values over the five runs listed in table III. 

The error in the holding time has two sources. First, during pressurization, the N 2 liquifies and diffusion 
begins before the final pressure of 1700 psig is reached. Second, after the final pressure is reached and high- 
pressure valve 3410N is closed, the pressure drops back as GN 2 continues to condense. The additional time for 
readjustment of the pressure adds to the holding time. The uncertainty in the holding time is estimated to be 
±0.5 min. Therefore the error in tf t is estimated as 


6th 

h 


0.5 

15 


x 100 = 3.3% 


(28) 


The error in flow rate R consists of three parts, as determined from equation (19): 


SR 

= 

6 R 0 

1 

— 

6 Pi 

1 

“1“ — 

6 T\ 

R 


R 0 

2 

Pi 

2 

Ti 


(29) 


The error in the reference flow rate is taken from the calibration sheet as 43 cm 3 /min. The errors in temperature 
and pressure at the flowmeter are due to fluctuations in ambient conditions, estimated to be ±5 out of 760 torr 
and ±3 out of 297 K. Then the error in R becomes 


SR 

~R 


f 43 1 5 1 3 \ 

\30 000 + 2 X 760 + 2 X 2977 


x 100 = 0.98% 


(30) 


The error in vi stems from errors in the handbook values of u/jy and vix and from the fact that xjvo 
is different for each of the runs (see eq. (24)). Since the former are considered accurate beyond the third 
significant figure, the error is attributed solely to variations in xjyo and is taken to be the standard deviation 
of the values of 1 listed in table III: 


6vl 

V L 


100 = 0.93% 


(31) 


The characteristic evaporation time T 2 is subject to two kinds of error: systematic errors, due to uncertainties 
in the O 2 analyzer reading and digitization error in data acquisition; and random error, due to deviations in 
system behavior from that assumed in the model. The latter are believed to outweigh the former by far. The 
error is taken to be the standard deviation of the values of T 2 listed in table III, which probably incorporates 
both kinds of error: 

&T2 
T2 


0.0708 

0.848 


x 100 = 8.35% 


(32) 


The surface area A of the liquid is known very precisely. The major source of error lies in the uncertainty 
in the 45° slant, which was determined to be ±0.025°. Thus 


SA 

A 


= (tan 45°) x 0.025° x x 100 = 0.44% 


( 33 ) 


Inserting the errors computed in equations (28)-(33) into equation (27) yields the total error in D: 


Pressure and Time Dependence of the LN2 Penetration 

Table IV summarizes the experimental results for Runs 6 and 7, conducted to determine the pressure 
dependence of the diffusivity, and Run 8, to determine the dependence of N 2 penetration upon holding time. 


Table IV. Summary of Results for Runs Used to Test the 
Pressure and Time Dependence of the N 2 Penetration 


Run 

N 2 pressure, 
psig 



zo/h 

v Li 

cm 3 /mol 

T 2i 

min 

6 

1000 

15 


0.9248 

29.327 

0.63 

7 


15 


1.0556 

28.263 

0.73 

8 


1 

0.170 

0.6743 

28.601 

0.19 


From the data given in the table, the diffusivities corresponding to N 2 pressures of 1000 and 500 psig are 
evaluated to be D = 4.9 x 10 -5 and 6.1 x 10~ 5 cm 2 /sec, respectively. Although these values are lower than 
those measured under an N 2 pressure of 1700 psig, the difference falls within the experimental error and is not 
considered significant. 

To investigate the time dependence, let the quantities associated with a holding time of 1 min (Run 8) be 
designated with a prime. Then, according to equations (7) and (9) 


II 

^1 




3.87 


(35) 


However, the measured ratio, according to the entries given in tables II and IV, is 


II 

r 2 


0.848 

0.19 


4.46 


(36) 


The difference between theory and experiment is 13%, well within the experimental error. This further confirms 
diffusion as the dominant mass transport mechanism. 

Impact on 8’ HTT Operations 

Now that the LN 2 -LOX diffusivity at 90.2 K has been established, it is possible to determine the LN 2 
contamination profile in the spherical 8000-gal LOX storage vessel of the modified 8’ HTT. Suppose it is 
specified that the LOX transferred from the storage vessel cannot have an N 2 contamination exceeding 1%. 
Then equation (1) can be used to predict how far below the LOX top surface the LN 2 concentration will exceed 
this limit for a given holding time: 

2 = 2(D t/ ! ) 1 / 2 erf _1 (l — 2x n) 

Using the value of diffusivity found earlier (see eq. (28)), and taking tj t = 15 min as a typical holding time, one 
finds 

2 = 2 x (8.6 x 10 -5 x 15 x 60) 1/2 x erf _1 (l - 2 x 0.01) = 0.9 cm 

Thus a layer at the top, about 1-cm thick, exceeds the contamination limit. If the top surface of the 
LOX happens to be at the equator of the storage vessel (worst case), the corresponding loss of LOX is 
430 000 cm 3 = 113 gal. The loss is insignificant compared with the total capacity of the storage vessel. 

Because the presence of any quantity of N 2 , no matter how small, lowers the boiling point of LOX, the 
contaminated LOX in the storage vessel is superheated after blowdown and will continue to evaporate, even 
without heat influx, as long as the pressure remains at 1 atm. Merely maintaining the pressure at 1 atm is 
an effective method of purging a contaminated layer from the main body of the LOX. After an acceptable 
contamination level is reached, e.g., < 1 percent, then it is essential to maintain a slight overpressure in the 
storage vessel to inhibit further evaporation. 
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Comparison With Prior Measurements 

Prior measurements of the LN 2 -LOX diffusivity were made by two Soviet investigators, Bahrov (ref. 6) at 
67.8 K and Boiko (ref. 7) at 77.8 K. Both used a version of the “gas phase” method, whereby gaseous N 2 
is added to O 2 vapor, which is in equilibrium with the liquid. The diffusivity is estimated from the rate of 
pressure change resulting from the diffusion of N 2 into the LOX. These and the present results are summarized 
in table V. 


Table V. Summary of Measurements of LN 2 -LOX Binary Diffusion 


Temperature, K 

— «» . . 

D, cm /sec 

Method 

Reference 

67.8 

1.07 X l0“ 5 ± 5% 

Gas phase 

6 

77.8 

2.8 x 10“ 5 ± 5% 

Gas phase 

7 

90.2 

8.6 x 10“ 5 ± 25% 

Differential flash vaporization 

This work 


The three sets of data listed in table V show a good fit to an expression of the Arrhenius form: 


D = Dq e 


AH 

Wl 


(37) 


where Tl is the temperature of the liquid. A plot of In D versus 1/X/, is a straight line, shown in figure 10 
together with the experimental points. A pressure correction is appropriate for the present measurement 
because the two other measurements were made at low pressures. However, such a correction is difficult to 
estimate because, in contrast to gases and solids, no theory of diffusion in liquids has proved successful in fitting 
the data well. A rough estimate based on a theory of Eyring (ref. 8), in which the diffusivity varies inversely 
with the viscosity, increases the present value of diffusivity by about 15%. The pressure dependence of the 
viscosity is taken from reference 4. Nevertheless, in view of the theoretical uncertainty and the inconclusiveness 
of the pressure dependence indicated by the present data, a pressure correction to the diffusivity plotted in 
figure 10 has not been applied. 



1/Temperature, K 1 

Figure 10. Arrhenius plot of the binary diffusivity: experiment (circles) and theoretical best fit (solid line). 
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A linear regression analysis yields the values for the preexponential factor Dq and activation enthalpy A H 
shown in table VI. The correlation coefficient is 0.99780, indicating an excellent fit. Table VI also shows values 
for two other cryogenic liquid systems (refs. 9 and 10), in which the diffusivity of a dilute radioactive tracer in 
LN 2 was determined. 

Table VI. Arrhenius Parameters for Binary Diffusion in Selected Liquid Systems 


Diffusing 

molecule 

Solvent 

Dq, cm 2 /sec 

AH, kcal/mol 

Method 

Reference 



2.38 x 10- ;i 

0.67 

Radioactive tracer 

9 

Ar 

n 2 

1.12 x 10“ 3 

0.59 

Radioactive tracer 

10 

n 2 

o 2 

4.52 x 10 -2 

1.08 

Differential flash vaporization 

This work 


The Arrhenius relationship expressed in equation (37) is characteristic of a thermally activated process. The 
associated physical model regards the liquid as having a short-range, quasi-crystalline structure with periodic 
potential wells. A diffusing molecule trapped in a well makes a certain number of attempts per second to jump 
out. The preexponential factor Dq is proportional to the attempt frequency, and the activation enthalpy AH 
is approximately equal to the height of the potential barrier. Inspection of the values of Dq listed in the table 
reveals that the attempt frequency for the N 2 -O 2 system is an order of magnitude greater than that for Ar or N 2 
in LN 2 , and the activation enthalpy is 1.5-2 times greater. The large difference in attempt frequency is difficult 
to explain theoretically and must be attributed to experimental error. The difference in activation enthalpy can 
be explained by a possible difference in mass transport mechanisms: Transport of a tracer molecule involves a 
single jump over a barrier, but the binary N 2 -O 2 transport could proceed by the double rotation mechanism 
of reference 11 and thus involve a double jump. 

Conclusions 

The differential flash vaporization technique has been adapted to determine the penetration of pressurized 
GN 2 into LOX, from which the binary diffusivity was evaluated. The mean penetration depth of five 
experimental runs, each at a GN 2 pressure of 1700 psig for a holding time of 15 min, was 0.9 cm. The 
corresponding binary diffusivity was found to be D = 8.6 x 10~ 5 cm 2 /sec ±25% at a liquid temperature 
of 90.2 K. Tests to determine the pressure dependence of the diffusivity were inconclusive. Another test 
confirmed, within experimental error, the expected square-root dependence of the penetration depth on time, 
thus providing additional support for diffusion as the dominant mass transport process. 

The measurement of the diffusivity obtained from this experiment, together with two prior measurements 
at lower temperatures, revealed an excellent fit to the Arrhenius equation, yielding a preexponential factor 
.Do = 4.52 x 10 -2 cm 2 /sec and an activation enthalpy AH = 1.08 kcal/mol. 

The contamination of LOX by pressurized GN 2 will have but minimal impact upon operations of the 
modified Langley 8-Foot High- Temperature Tunnel. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
February 14, 1989 
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Symbols Introduced in the Appendices 


A R correction to the evaporation rate, mol/min 

s state variable of the Laplace transform 

U flow rate of vapor leaving the ullage volume, mol/min 

X(s) Laplace transform of xj^(t) 


= t 
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Appendix A 

Derivation of the Time History of the N2 Molar Content 

Applying the Laplace transform to equation (10), one obtains 


sX{s) - Xm + ^ (S> = 2^7 ' 2^ £ { erf (j + 4 ) } 

where X(s) is C{xj^(t)}. Equation (Al) can be solved for X(s): 

X N0 , 1 C{eri{z 0 /h + t/r 2 )} 


X(s) = 


s + 1/t\ 2t\s(s + 1 /rj) 


2ri(s + l In) 


(Al) 


(A2) 


The inverse transform of the first two terms on the right are well known, but the difficulty is to find that of 
the third term. One proceeds by starting with the transform pair 5.12-4 in reference 12: 


erf ^ if -I- — erf (a) | 


(. s — a)~ l e s erfc(s) = C | e a ( t+ “) 

Replace s by s/a and apply the scaling law: 

- a) e ( s / a ) 2 er f c = C j ae a ^ at+a ' > erf — erf (a) | 

Now replace s by s + aa and apply the translation law: 


ae [(s+aa)/a) erfc (£±pa) 
(s — aa + aa) 


= £ {““ V “ ,e “‘H(f-)-erf W ]} 


which can be simplified to 


- = £ j err (f + a) - erf(a)| 


Let a = zo/h and a = 2/t 2 and substitute into equation (A6): 

e _ ( z o/^) 2 e( sr 2/2+zo/^) 2 er f c ( tS7 - 2 /2 + zq/K) erf(za/h) 




)} 


Substitute equation (A7) into (A2) and expand the last term into partial fractions: 


X(s) = 


X N 0 


+ 


1 — er{(zo/h) 


(ST 2 


( l - 1 ^ 

V 2 

h h) 

s + I/tiJ 


s + 1/ri 2ris(s + 1/ri) 2 

To evaluate the last term replace s by s + (3 and apply the translation law to equation (A6): 

erfcK.+w«w;i = c i. j erf ^ + _ erf(al) j j 

Let /3 = A and a 1 = ^ ^ and substitute into equation (A9): 


^/2+^o/fe) 2 erfc (g r2 / 2 + ZQ /h) 

s + 1/ri 


(A3) 


(A4) 


(A5) 


(A6) 


(A7) 


(A8) 


(A9) 


= cfe^e^ °/ ft -V 2 n ) 2 Lf (- + ^ - erf ( ^ - — ')] } 

l Vr 2 h 2 tJ \h 2 t\ )\ j 


(A10) 
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Upon applying equations (A7) and (A10) to evaluate the last terms, one finds the time history xjy(i) b 
the inverse transform of equation (A8) and substituting equation (11) for x'.yo . The solution is 


x N {t) = 


1 


1 - erf ( - + ^ 
T2 h 


+ iexp 


L + _ II. 

ti hri 4 7-j 


f , t zq r 2 

er M b T 7 , — 

T2 h 2 t\ 


— erf 


£0 

h 


taking 



(All) 
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Appendix B 

A Method of Correcting for a Slow Pressure Buildup in the Ullage Volume 


A correction to the evaporation rate due to a slow pressure buildup in the ullage volume is derived in the 
following. From the continuity equation 

—A* 


where U is the flow rate of vapor leaving the ullage volume (1.23 mole/min), which was held constant by 
adjustment of valve A. The factor \/2 accounts for the slant of the run tank . From equation (15), parameters 
V and T are assumed to remain constant, but P builds up very slowly with time in the ullage volume. Therefore, 
equation (Bl), combined with equation (15), can be expressed in the following form: 


R = U + 


1 d f PV \ _ 
dt \jwf ) ~ 


v/2 


U + 


V dP 
s/2R Q T~dt 


(B2) 


Substitute the values V = 8000 cm 3 , R g = 82.1 cm 3 -atm/mol-K, and T = 90.2 K into equation (B2) to obtain 
a correction to the evaporation rate: 

dP 

AR = R-U = 0.76388— (B3) 


dt 


The time derivative of P follows from equation (16): 

dP 


— = 6 2 + 63(61 - l)e~ M 


(B4) 


where the pressure has been converted to atmospheres. Upon substituting the values of b\, 63 (see table Bl), 

and t = (including the offset time) into equation (B4), one can evaluate A R. 


Table Bl. Summary of Results for the Correction of T 2 Due to Pressure Buildup 


Run 

61 

62 

63 

min 

ft 

min 1 

(*)„• 

atm/min 

^ x 100% 

1 

3.1411 



0.92475 

-0.03547 



2 

6.86 



1.33575 

-0.02120 


liSscSM* 

3 

3.255 



1.33575 

-0.01631 



4 

3.959 

0.1099 


1.13025 

-0.01916 

0.430 

-5.49 

5 

3.973 

0.18344 

Sftu9 

0.92475 

-0.03609 

0.713 

-17.3 


Next, the modified transport equation due to pressure buildup is determined by again substituting equations 
(2), (7), (8), and (9) into equation (6) and changing the evaporation rate R on the left-hand side to the flow 
rate of vapor U : 


dxjy 

dt 


1 £ 

-I x N - — 

n 2ri 


1 — erf 


OK)] 


(B5) 


where £ = R/U and t\ is redefined as ny/U. 
Differentiation of equation (B5) leads to 


^ X N 1 dxjq _ £ ,-(zn Ih+tlroS 2 

rif 2 T\ dt \pKT\T 2 


(B6) 


At the inflection point M , 



= 0. Thus, equation (B6) becomes 


T2e (tM/r2+zo/h) 2 _ h 

y/n(dxff/dt) m 


(B7) 
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By differentiating equation (B7), one obtains the correction d,T 2 in the following form: 

1 e -{iMl T 2+^/h) 2 d c 

dr 2 = jj — r — y~ z — 77 (B8) 

^(^)m I 1 - ^ ft + ?)] 

where d£ = AR/U. The value of dr 2 is found by substituting the parameter values obtained from tables III 
and B1 into equation (B8). The percentage corrections to r 2 due to pressure buildup for the 1700 psig/15 min 
runs are listed in the final column of table Bl. The relative correction to T2, however, is still smaller than the 
error of T2 (except for Run 5); thus it is incorporated in the systematic and random errors of 72 . Hence, the 
pressure buildup in the ullage volume did not considerably affect the experiment. 
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